Cosmic String Loops Collapsing to Black Holes 

"O ! R.N. Hansen* M. Christensen"^ and A.L. Larsen^ 

O 

^ ■ February 7, 2008 

>.■ 

^ : 

Physics Department, University of Odense, 
VO I Campusvej 55, 5230 Odense M, Denmark 

(N 
> 

oo 

o 

CN ■ Abstract 

o 

0\ 



We examine the question of collapse of Turok's two-parameter fam- 
ily of cosmic strings. We first perform a classification of the strings 
O^' according to the specific time(s) at which the minimal string size is 

reached during one period. We then obtain an exact analytical expres- 
sion for the probability of collapse to black holes for the Turok strings. 
Our result has the same general behavior as previously obtained in the 
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H ' literature but we find, in addition, a numerical prefactor that changes 

the result by three orders of magnitude. Finally we show that our 
careful computation of the prefactor helps to understand the discrep- 
ancy between previously obtained results and, in particular, that for 
"large" values of G/i, there may not even be a discrepancy. We also 
give a simple physical argument that can immediately rule out some 
of the previously obtained results. 
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1 Introduction 

According to the hoop-conjecture 0, a cosmic string that contracts to a size 
smaller than its Schwarzschild radius will collapse and form a black hole. This 
process is of particular interest and importance in connection with primordial 
black holes, Hawking radiation, high energy cosmic gamma bursts etc. 

In a network of cosmic strings, only a very small fraction, /, of strings 
are expected to collapse to black holes. Many attempts have been made to 
obtain a value for / (see for instance 0-||T^), but the results deviate wildly. 



Interestingly enough, only in one of the pioneering papers on the subject, 
the one by Polnarev and Zembowicz |^ , is the derivation of / based on exact 
analytical expressions for the cosmic strings involved in the process. In all 
other discussions the derivation of / is based on linearized expressions for 
the string configurations (for instance [0, 3), rather general arguments and 
estimates (for instance [^ ^]), observational data concerning high energy 
cosmic rays (for instance 0-0) or the derivation is purely numerical (for 
instance [|10|, |ll|). And even in Ref. 0, the final computation of / is actually 
numerical, although it could in fact have been performed analytically. 

In the present paper we consider the question of collapse for the analytical 
two-parameter family of strings introduced by Turok |TB|. This is the same 
family of strings that was considered in p[ . However, simple explicit examples 
show that the results obtained in [Q are correct only in part of the parameter- 
space. We will show that some essential points were missed in Ref. 0, and 
that the results obtained there are not completely correct. 

First of all, we make a classification of the string configurations according 
to their general behavior during one period of oscillation. Together with 
some simple explicit examples, this analysis reveals that the corresponding 
result obtained in 0] is in fact only correct in approximately half of the two- 
dimensional parameter-space. Secondly, we then derive the exact analytical 
expression for the probability / of string collapse to black holes. Our result 
for / agrees partly with that of Ref. in the sense that / oc (G/x)^/^, 
where n is the string tension and G is Newtons constant. However, we find 
a numerical prefactor in the relation, expressed in terms of Euler's gamma- 
function, which is of the order 2000. For comparison, in Ref. the prefactor 
was found (numerically) to be of the order 1. Finally we show that our careful 
computation of the prefactor helps to understand the discrepancy between 
previously obtained results and, in particular, that for "large" values of G/i, 
there may not even be a discrepancy. We also give a simple physical argument 



that can immediately rule out some of the previously obtained results. 

The paper is organized as follows. In Section ^, we introduce the string 



configurations of Turok |T5[ and discuss the results obtained by Polnarev and 
Zembowicz p[. In Section |^, we give the precise definition of the minimal 
string size R, and we make a complete classification of the two-dimensional 
parameter-space. The classification is in terms of the specific time(s) at 
which the value R comes out for a particular string configuration during 
one period of "oscillation". In Section ^, we derive an exact analytical 
expression for the probability of string collapse to black holes and we give 
the approximate result in the physically realistic limit Gfi << 1. Finally 
in Section ||, we argue that our careful computation of the prefactor helps 
to understand the discrepancy between previously obtained results, and we 
give our concluding remarks. Some of the details of the computations used 
in sections |^ and ^ are presented in the Appendix. 



2 Two- Parameter Turok-Strings 

The string equations of motion in fiat Minkowski space take the form: 

X^-X"^ = (2.1) 

supplemented by the constraints: 

V.uX'^X''^ = T],, [X'^X'' + X'^^X''^) = (2.2) 

It is convenient to take X^ = At, where A is a constant with dimension of 
length. Then Eqs. (|2.1| )-( p^ become: 

i = X" 

X-X' = (2.3) 



X2 + X 



A' 



A two-parameter solution to Eqs. ( p.3|) was first introduced by Turok |T5 
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cos3(a - r) + (1 - 2/5) cos(or + r) 



2^Ja{l -a)cos{a - r) + 2^ (3{l - f3)cos{a + r) 
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(2.4) 



where a G [0; 1], /9 G [0; 1]. This family of solutions generalizes the solutions 
considered in Ref. [p!6| . 

The total mass-energy of the Turok-string is: 

Energy =- f d^X T\ (2.5) 

where: 

T^, = -L^ = /i I drda rj^xrj.s {x^X' - X'^X'') 5^ (X - X(r, a)) 

(2.6) 
It follows that: 

Energy = 2TTAfi (2.7) 

which (by construction) is independent of the parameters (a,/?). Similarly, 
one finds for the momenta: 

'^3x^0 = (2.8) 

In fact, the string center of mass is located at Xcm = (0, 0, 0) at all times, 
and the strings are symmetric under reflection in origo. The Schwarzschild 
radius corresponding to the energy (^TTf) is 



Rs = AnAGfi (2.9) 

For realistic cosmic strings, the dimensionless parameter Gfi is quite small 
|13, m, say Gfi ~ 10"^. It is then clear that a Turok-string (|]D will 
typically be far outside its Schwarzschild radius (as follows since typically 
X ^ A, Y r^ A, Z r^ A). However, for certain particular values {a, (3), 
a string might at some instant during its evolution be completely inside its 
Schwarzschild radius. Such a string will collapse and eventually form a black 
hole, according to the socalled hoop-conjecture |I|. To determine whether 
this happens or not, we must first find the minimal 3-sphere, that completely 
encloses the string, as a function of time. After minimization over time, we 
then get the minimal 3-sphere that can ever enclose the whole string. Let 
the radius of this sphere be R (it will be defined more stringently in the next 
section). Then the condition for collapse is: 

R<Rs (2.10) 



That is to say, 

^<47rG/i (2.11) 

In a pioneering paper, Polnarev and Zembowicz considered, among 
other things, the question of collapse of the Turok-strings ( ^.41) . In connection 
with the minimal string size R they found: 



The strings have their minimal size R at 



vr 
2 



(2.12) 



For generic parameters (a,/3): 



The result ( p.l2| ) would be expected for a string experiencing a monopole- 
like oscillation, i.e. starting from maximal size at r = 0, then contracting 
isotropically to minimal size at r = 7r/2, and then re-expanding to its original 
maximal size at r = tt. As for the result ( |2.13D it was simply stated 
without proof or derivation. 



In Sections 3,4 we shall show that ( |2.12| ) and ( |2.13D are not completely 



correct. In fact, they are correct in part of the two-dimensional parameter- 
space but incorrect in other parts. A simple example showing that (|2.13|) 



cannot be correct, is provided by the case {a, (3) = (3/4, 1/4). In that case 
the result (p. 13 ) would actually give R = 0, which would mean that the 



whole string had collapsed to a point at some instant. This would imply 
that at some instant X = Y = Z = 0, but that is certainly impossible for 
(a,/3) = (3/4,1/4). The problem is that ( ^.13| ) gives the minimal distance 



from origo to the string, but it is the maximal distance which is relevant for 
the minimal 3-sphere. 



3 Minimal String Size and Classification 

For a given pair of parameters {a, (3), we define the minimal string size R^ 
as the (square of the) radius of the minimal 3-sphere that can ever enclose 



the string completely. More precisely: 

(3.1) 



2 _ Minimum 

r G [0; vr] 



Maximum / 2 

a e [0; -K 



[R^{T,a 



where: 

R^ (r, a) = X\t, a) + ^'(r, a) + ^^(r, a) (3.2) 

as obtained using Eq. (|2.4|) . Thus, for a fixed time r, we first compute 
the maximal distance from origo (i.e. the string center of mass) to the 
string. This gives the minimal 3-sphere that encloses the string at that 
particular instant. We then minimize this maximal distance over all times. 
This gives altogether the minimal string size R^. And this is obviously the 
quantity that must be compared with the (square of the) Schwarzschild 
radius -R|. Notice that we need only maximize over a G [0; vr] and minimize 
over r G [0; tt] in Eq. (|3.1|) . This is due to the reflection symmetry and time 
periodicity of the Turok-strings (p.4|). 



We now outline the computation of R^. The details can be found in the 
Appendix, and some analytical results are given also in Section ^. 
We first solve the equation: 

^^ = (3.3) 

for fixed time r. This leads to a quartic equation in cos(2cr): 

cos^(2a) + a cos^(2or) + b cos2(2ct) + c cos(2a) + d = (3.4) 

where the coefficients (a, b, c, d) depend on time r as well as on the parameters 
{a,f3). The explicit expressions are given in the Appendix. The solutions 
to equation ( |3^ ) are explicitly known, leading to cr = a^r) for given values 
of (a,/?). By insertion of these solutions a = cr(r) into -R^(r, a), it is then 
straightforward to obtain the maximal distance in the square bracket of Eq. 
( |3.1| ). This is now a function of r, which finally has to be minimized over 
r G [0; vr]. For more details, see the Appendix and Section ^. 

The result is shown in Fig. 1, that is, the minimal string size R^ as a 
function of [a, (3). Not surprisingly, the minimal string size is generally of 
the order 1 (in units of A; see Eq. ( p.4| )). The exception is in the vicinity of 
{a, (3) = (0,0), where the minimal string size is close to zero. This was also 

6 



to be expected since the vicinity of (0, 0) describes the near-spherical strings, 
and those are the only ones expected to have a chance to collapse, due to 
their relatively small angular momentum. These strings will be considered 
in detail in Section ^ 

However, the computation of the minimal string size R^ gives some 
more information about the Turok-strings, namely the time(s) r at which 
the strings are minimal. This gives rise to a precise classification of the 
Turok-strings, and a subsequent subdivision into 3 different families (see 
Fig. 2): 

I. These strings have their minimal size at r = 7r/2. That is, starting from 
their original size at r = 0, they generally contract to their minimal size at 
T = n/2, and then generally expand back to their original size at r = vr. 

II. These strings start from their minimal size at r = 0. Then they generally 
expand towards their maximal size and then recontract towards their 
minimal size at r = vr. 

III. These strings have their minimal size at two values of r symmetrically 
around 7r/2. That is, they first generally contract and reach the minimal size 
at some Tq G [0;7r/2]. Then they expand for a while, and then recontract 
and reach the minimal size again at r = vr — tq. Then they expand again 
towards the original size at r = vr. In this family of strings, the value of Tq 
depends on (a,/3). 

It must be stressed that the strings in most cases do not expand or 
contract isotropically. They typically expand in some directions while con- 
tracting in other directions. This is why we use the expressions "generally 
expand" and "generally contract" , which refer to the minimal string size as 
a function of time (the radius of the minimal 3-sphere enclosing the string, 
as a function of time) . 

Notice that besides the three above-mentioned families of strings, there 
are a number of degenerate cases at the different boundaries. For instance, 
at the boundary a = 1, the strings have their minimal size at r = 0, 
r = 7r/3, r = 27r/3 and r = vr. On the other hand, at the boundary between 
regions I and II, the strings have their minimal size at r = 0, r = 7r/2 and 
T = 71. Notice also that the 3 points {a, (3) = (1,0), {a,j3) = (1,1) and 
[a, (3) = (0,1) correspond to rigidly rotating strings, thus they have their 



minimal (and maximal) size at gM times. 



Let us close this section with a comparison with the result (|2.12|) of Ref. 



We see that the result ( p.l2| ) is correct in the region I of parameter-space, 



but incorrect in regions II and III. 



4 Probability of String Collapse 

In this section we consider the question of collapse of the Turok-strings 
As already discussed in the previous section, the only strings with a chance to 
collapse are those corresponding to parameters (a,/3) located in the vicinity 
of (0,0). That is, we need only consider strings in the family I of Fig. 2. 
Using the results of the Appendix, it is then straightforward to show that 
the minimal string size R^, as defined in Eq. (|3.1|), is given by: 



where: 



R^ = Max(^Rl Rl) (4.1) 



^_4a2 
A^~ 9 
Rl ( r—. r nrr. ^X' fa 



(4.2) 



J=(V"(l-«)-\//5(l-/5)) +(3-^) (4.3) 

Notice that Eq. ( [4.3| ) is precisely the result ( p.l3| ) of Polnarev and Zembowicz 
p|. However, in Ref. 0, the other solution ([4.2|) was completely missed, and 
this is actually the relevant solution in Eq. ( |4.1| ) in approximately half of the 
parameter- space (a,/3). 

According to the hoop-conjecture (see for instance ||1||), the condition for 
collapse to a black hole is then given by Eq. ( p.ll| ), with R given by Eqs. 

(E3)-(iD: 



Max (rI, 



Ri) < {AiiAG^iY (4.4) 



which should be solved for (a,/?) as a function of G/i. This can be easily 
done analytically. The result is shown in Fig. 3: the part of parameter-space 
fulfilling inequality ([4.4|) is bounded by the a-axis, the /5-axis, the straight 
line a = 3Rs/2A and the two curves: 



^ , , -16a^-12a^ + 27a-9(2a-3)Rl/A^±6VD ,^ , 
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where: 

D = —a (1 — a) 
Notice also that, 



A J A , 



/3+(0) 



-^1 

A2 



ia- - 6« + 9 ('^ + l) ^ 



(4.6) 



and 



P- («o) = 



for ao = — I 1 



l^_32Rl 



9 A^ 



(4.7) 



The probabihty for collapse into black holes is then given by the fraction /: 



JR 



R<Rs 



dad(3 



"0 rP+{oi) 

da I d(3 + 

Jo Jao 

3-Rg 
2A 



^ . rM») 



da I d(3 

/3_(a) 



/9+ (a) da + [/?+ (a) — /?_ (a)] da 



(4.8) 



This equation represents the exact analytical result for the probability of 
collapse of the Turok-strings ( |2.4| ), for a given value of Rs/A = iTiGfi. The 
integrals in ( [4.8| ) are of hyper-elliptic type |T9[, and not very enlightening 
in the general case. However, using that typically G/i -C 1 (see ||l^ HH), a 
simple approximation is obtained by keeping only the leading order terms: 



, 12^6,, „ 5 /•! t'^dt ,^(,r^ a 

f = ^— (47rG/i)2 I -^== + O ({G^iY 



5r2 



^(G/i)t +0((G/i)i) 



(4.9) 



The result ( [4.9| ) is a very good approximation for G^ < 10 ^, thus for any 
"realistic" cosmic strings we conclude: 



/^2-103-(G/i)^ 



(4.10) 



which is our final result of this section. 

It should be stressed that we have been using the simplest and most 
naive version of the hoop-conjecture: namely, we did not take into account 



the angular momentum of the strings. However, numerical studies |]TT| of 



other families of strings showed that inclusion of the angular momentum 
only leads to minor changes in the final result, so we expect the same will 
happen here. It should also be stressed that we have neglected a number of 
other physical effects that might change the probability of collapse. These 
include the gravitational field of the string and gravitational radiation. 

Finally, as in all other discussions of the probability of string collapse, we 
are faced with the problem that we do not know the measure of integration 
in parameter-space. Thus using another measure in Eq. ( [4.8|) would 
generally give a different result (see also Ref. 0). 



5 Conclusion 

In this paper we examined, using purely analytical methods, the question of 
collapse of Turok's two-parameter family of cosmic strings |1^. We made 
a complete classification of the strings according to the specific time(s) the 
minimal string size is reached during one period. This revealed that the pre- 
viously obtained results |0] were only correct in part of the two-dimensional 
parameter-space. 

We then obtained an exact analytical expression for the probability / of 
collapse into black holes for the Turok strings, which partly agrees with that 
of Ref. [0 in the sense that / oc (G/x)^/^. However, we showed that there is 
a large numerical prefactor in the relation. This factor is of the order 2000, 
and not, as previously stated [0, of the order 1. 

One might say that it is perhaps not so important whether the prefac- 
tor is 1 or 2000 since the exponent will more or less kill the probability of 
collapse anyway. This may very well be true for "small" values of Gfi (say, 
G/i ~ 10^^), but for "large" values of Gfi (say, Gfi ~ 10^^) the situation 
is completely different. In fact, we shall now argue why it is so important 
to carefully compute the prefactor: We find that when using our result, a 
clear picture is beginning to emerge. Different computations based on differ- 
ent families of strings (not surprisingly) produce slightly different exponents, 
but they also produce completely different prefactors. Importantly, the two 
go in the same direction: An exponent larger by 1 is followed by a prefactor 
larger by a factor 100 (roughly speaking). For instance, the one-parameter 
family of Kibble and Turok [^] gives / ~ 167r^(Gyu)^, our computation for 



the two-parameter family gives / ~ 2000(G/i)^/^ and the Caldwell-Casper 
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computation [|n[] gives / ~ 10^ (Gfi)^. We find it extremely interesting that 
for Gfi ~ 10~^ (which is the range where the Caldwell-Casper computation 
is vahd), the three computations basically agree giving / fa 10~^ — 10~^. We 
therefore find that our careful computation of the prefactor is very important 
and has helped to understand the discrepancy between previously obtained 
results. And in particular, for "large" values of Gfi, there may not even be a 
discrepancy since the different exponents and prefactors of the different com- 
putations actually produce more or less the same number for the probability 
of collapse. For "small" values of G/i, the picture is unfortunately not so 
clear, and more detailed investigations seem necessary. 

It is actually possible to give a physical argument showing that the results 
of different computations (if they are done correctly) should merge for "large" 
values of Gfi. Consider the rigidly rotating straight string, corresponding to 
a = and P = 1. It is easy to show that it is precisely inside its Schwarzschild 
radius for G/j, = (47r)~^. Now, it is well-known that the rigidly rotating string 
has the maximal angular momentum per energy (it is exactly on the leading 
Regge trajectory), and therefore is expected to be the most difficult string 
to collapse into a black hole. Therefore, for Gn ~ (47r)~^ all other strings 
have already collapsed, and we should expect / ~ 1. This is indeed the case 
for the three above mentioned computations, while it does not hold for the 
computation of Ref.|0. More generally, the condition /((47r)^^) ^ 1 can 
be considered as a physical boundary condition, and therefore can be used 
to immediately rule out some of the previously obtained results for /; for 
instance those of [|| and p. 

As a possible continuation of our work, it would be very interesting to 
consider more general multi-parameter families of strings, to see how general 
our result for / actually is. Such families of strings have been constructed 
and considered for instance in |2y, |^ , and more general ones can be obtained 
along the lines of [^ . 



Unfortunately, there are also still some open questions, as we discussed 
at the end of Section ^. The main problem seems to be that we still do not 
know exactly what is the measure of integration in parameter-space H] . 
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Appendix: The Minimal String Size 

In this appendix, we give some details of the results presented in the sections 
I and |. 

The distance from origo to the string, as a function of r and a, is con- 
veniently written as a polynomial in sin (2cr) and cos (2cr) . From (|3.2| ) and 
(121): 



A^ 



Aa/3 



R" (r, a) = _{Co + — ^ [Ci sin {2a) + C2 cos {2a)] 



Aaf3 



cos (2r) cos (2a) + sin (2r) sin {2a)] cos {2a) } (A.l) 



where: 



Co = 


- 2-^ 


^ + 2[2^a{l-a)^P{l-P) 




+ (l-a)(l-/3) + ^]cos(2r) 


Ci = 


3 


"^(l-/3)sin(4r) + ^(l-a)sin(2r)' 


C2 -- 


3 

Aa(3 


^(l-/?)cos(4r) + ^(l-a)cos(2r) 




+ 


2[2^a{l-a)^P{l-P)-P{l-a)J 



(A.2) 



Then it is straightforward to show that the condition (|3.3|) leads to 

x^ + ax^ + bx'^ + cx + d = (A. 3) 
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where x = cos {2a) and: 



a 
b 



-Ci sin (2r) - C2 cos (2r) 
Cf + CI ^ 



(A.4) 



Ci 



sin (2r) + C2 cos (2r) 



sin^ (2r) - Ci 



The solutions to Eq. (|A.3|) can be written down in closed form. Define 



y 
z 



6^ - Sac + I2d 

2b^ - 9abc + 27c^ + 27a^d - 72bd 

1 



(A.5) 



23^" Z 

w = —— + 



3Z 



3 2^ 



Then the 4 solutions are: 



2:1,2 



2:3,4 



a 

'4 



1 a2 
2VT 



26 

y 



Wt 



'\ 



a^ 
~2 



4b 



W 



-a^ + Aab - 8c) 



2b 



+ w 



a 1 a2 2b ,,, 1 

h -\ h W=F - 

42V 4 3 ^2 



\ 



!-- 



(-a3 + Aab - 8c) 

Vf-f +"^ 



(A.6) 



which give a = cr (t) for given values of (a,/3). These solutions are 
inserted into ( |A.1|) and then the result of the square bracket in Eq. ( p.l|) is 
determined. Finally this function of r must be minimized. 



As an example, consider strings in the region I of parameter-space; see 
Fig. ^. Using the above formulas, one finds: 



2:1,2 
2:3,4 



_C2 

2 



(A.7) 
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with 



Co 



2a/3 



(f + /5)(«-/3)-(V«(l-«)-\//5(l-/5))" 



(A.8) 



Insertion into Eq. ( [A. 1| ) then leads directly to the result for the minimal 
string size in the region I, 



i?2 = Max(i?^, Rl) 



where 



AA'a 



2^,2 



-"^-l - 


c 








^ ? 


\ 21 


i?2 - 


= A' 


(v/a (1 - 


-")- 


-\//5(l- 




-.)] 



(A.9) 



(A.IO) 



c.f. Eqs. (PD-Q 
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CAPTIONS FOR FIGURES 



Figure 1: The radius of the minimal 3-sphere completely enclosing a string 
with parameters (q;,/3) plotted for all parameter-space. Notice that R^ is 
close to zero only near {a, (3) = (0,0), that is for the near-spherical string 
configurations. 



Figure 2: The considered strings fall into three families. The ones that reach 
their minimal size i? at r = 7r/2 (I), at r = 0, tt (II) and at r = tq, vr — tq 

forroe]0,7r/2[(III). 



Figure 3: The region of parameter-space which contains the strings falling 
inside their own Schwarzschild radius is bounded hj a = 0, f3 = 0, a = 
3Rs/2A and the two curves l3±. Here Gfi = 10~^ is chosen in order to 
illustrate the general form of the region. As Gfi decreases, the region is 
relatively stretched out and becomes quite narrow. 
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